Feedback control of a solid-state qubit using high-fidelity projective measurement 
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We demonstrate feedback control of a superconducting transmon qubit using discrete, projective 
measurement and conditional coherent driving. Feedback realizes a fast and deterministic qubit 
reset to a target state with 2.4% error averaged over input superposition states, and cooling of 
the transmon from 16% spurious excitation to 3%. This closed-loop qubit control is necessary for 
measurement-based protocols such as quantum error correction and teleportation. 
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Many protocols in quantum information processing 
(QIP) require closing a feedback loop where coherent con- 
trol of qubits is conditioned on projective measurements 
in real time [1]. Important examples include quantum 
error correction and teleportation [2], so far achieved 
in ion-trap [3, 4] and photonic systems [5, 6]. During 
the last decade, the steady development of qubit read- 
out and universal gates needed in a quantum proces- 
sor [7] has made superconducting circuits [8] a leading 
solid-state QIP platform. However, the simple quantum 
algorithms [9] and teleportation-like protocol [10] so far 
demonstrated fall in the category of open-loop control. 
Measurement is performed as the final step, following a 
programmed sequence of applied gates. A comparable 
realization of closed-loop control has been precluded by 
stringent requirements on high measurement fidelity and 
short loop delay (latency). Until recently, the available 
qubit coherence times bottlenecked both achievable fi- 
delity and required speed. 

For feedback control of superconducting qubits, the 
development of circuit quantum electrodynamics [11, 12] 
with 3D cavities (3D cQED) [13] constitutes a watershed. 
The new order of magnitude in qubit coherence times 
(> 10 /is) allows boosting projective-readout fidelity up 
to 98% [14, 15] and realizing feedback with off-the-shelf 
electronics. Very recently, feedback based on continuous 
weak measurement has sustained Rabi oscillations of a 
transmon qubit indefinitely [16]. Previously, this type 
of feedback had only been used to generate and stabilize 
quantum states of photons [17], ions [18], and atoms [19]. 

In this Letter, we demonstrate feedback control of a 
superconducting transmon qubit based on discrete, pro- 
jective measurement. This dual type of feedback is the 
kind necessary for measurement-based QIP. As a first ap- 
plication, we demonstrate a feedback-based reset that is 
deterministic and fast compared to passive initialization. 
This feedback cools the transmon from a spurious steady- 
state excitation of 16% to 3% and resets qubit states 
with 2.4% error averaged over the Bloch sphere. These 
absolute errors are dominated by latency, in quantitative 



agreement with a model including transmon equilibration 
and readout errors. 

The experiment employs a transmon qubit inside 
an aluminum 3D cavity [13]. The qubit (w i/27r = 
4.889 GHz transition frequency) couples to the cav- 
ity fundamental mode (u> t /2tt = 6.546 GHz, coupling- 
limited linewidth k/27t = 550 kHz) with strength g/2ir = 
68 MHz. The high-fidelity, projective qubit readout 
forming the input to the feedback loop uses homodyne 
detection of the qubit-state dependent cavity transmis- 
sion (dispersive shift 2x/2n = -1.9 MHz [12]). A 400 ns 
measurement pulse at wrf = w r — X is applied to the 
cavity and the transmitted signal is then amplified by a 
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FIG. 1. (color online). Transmon equilibration to steady 
state. Time evolution of the ground-state population P| ) 
starting from states |0) and |1). Solid curves are the best 
fit (including data in Fig. S4) to Eq. 1, giving the inverse 
transition rates T^ 1 = 50 ± 2 fis, — 20 ± 2 fis, T^, 1 = 
324±32 /is, r^ 1 = 111±25 /is. From the steady-state solution, 
we extract residual excitations f|i) ]SS = 13.1 ± 0.8%, P|2),ss — 
2.4±0.4%. Inset: steady-state population distribution (bars). 
Markers correspond to a Boltzmann distribution with best-fit 
temperature 127 mK, significantly higher than the dilution 
refrigerator base temperature (15 mK). 
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Josephson parametric amplifier [20, 21] to enhance sen- 
sitivity [22]. The feedback controller is an ADwin Gold 
processor that samples the transmitted homodyne sig- 
nal, performs 1-bit digitization to interpret the projected 
qubit state, and conditionally triggers a tt pulse resonant 
with the transmon f> 1 transition. The 2.64 /its delay 
between start of the measurement and end of the ir pulse, 
set by processing time in the ADwin, is short compared 
to the qubit relaxation time T\ (see below) [23]. 

Our first application of feedback is qubit initializa- 
tion, also known as reset [7]. To be useful in QIP, re- 
set must be deterministic (as opposed to heralded or 
postselected [14, 15]) and fast compared to qubit coher- 
ence times. Obviously, the passive method, i.e., waiting 
several T\ times, does not meet the speed requirement. 
Moreover, it can suffer from residual steady-state qubit 
excitations [14-16, 24], whose source remains an active 
area of research. The drawbacks of passive initialization 
are evident for our qubit, whose ground-state population 
P|o) evolves from states |0) and |1) as shown in Fig. 1. 
The transmon dynamics are captured by a master equa- 
tion model for a three- level system [22, 25]: 

^io> \ / -r i0 r 01 o \ / p lQ) \ 
P\i) = r 10 -r i-r ai r 12 ^ . 

P\2) J V ^1 / V P\2) J 

(i) 

The best fit to the data gives the qubit relaxation time 
Ti = l/r 01 = 50 ± 2 fis and 15.5% residual total exci- 
tation. Previous approaches to accelerate qubit equili- 
bration include coupling to dissipative resonators [26] or 
two-level systems [27]. However, these are also suscepti- 
ble to spurious excitation, potentially inhibiting complete 
qubit relaxation. 

Feedback-based reset circumvents the equilibration 
problem by not relying on coupling to a dissipative 
medium. Rather, it works by projecting the qubit with 
a measurement (Mi, performed by the controller) and 
conditionally applying a ir pulse to drive the qubit to a 
targeted basis state (Fig. 2). A final measurement (M 2 ) 
determines the qubit state immediately afterwards, by 
averaging the second half (200 ns) of a transmitted mea- 
surement pulse sampled by a data acquisition card [22]. 
In both measurements, the result is digitized into levels 
H or L, associated with |0) and |1), respectively. The 
digitization threshold voltage I4h maximizes the readout 
fidelity at 99%. The ir pulse is conditioned on Mi = L to 
target |0) (scheme Fb ) or on Mi = H to target |1) (Fbi). 
To benchmark the reset protocol, we quantify its action 
on the purest states we can prepare. This step is accom- 
plished with a preliminary measurement M (initializing 
the qubit in |0) by postselection [14, 15]), followed by a 
calibrated pulse resonant on the transmon 4-> 1 transi- 
tion to prepare \6) = cos(0/2) |0) + sin(0/2) |1). We first 
consider the action of reset on the basis states (9 — 0, tt). 
The overlap of the M 2 histograms with the targeted re- 
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FIG. 2. (color online). Reset by measurement and feed- 
back, (a) Histograms of 300 000 shots of Mi, with (squares) 
and without (circles) inverting the qubit population with a n 
pulse. H and L denote the two possible outcomes of Mi, dig- 
itized with the threshold Vth, maximizing the contrast. Full 
(empty) dots indicate (no) postselection on Mo > V pa , chosen 
to prepare |0) with 0.7% error while keeping 81% of the data. 
This technique purifies the steady-state mixture (see Fig. 1), 
preparing the input states for the feedback sequences in (b) 
and (c). (b) Histograms of M2 after applying the feedback 
protocol Fbo, which triggers a tv pulse when Mi = L. Using 
this feedback, ~ 99% (92%) of counts fall in the H area for 
9 — (n), respectively, (c) Feedback with opposite logic Fbi 
preparing the excited state. In this case, ~ 98% (94%) of 
counts fall in the L area for 8 — (tt). 

gion (H for Fbo and L for Fbi) averages at 96%, indicat- 
ing the success of reset. Imperfections are more evident 
for 6 = tt and mainly due to equilibration of the trans- 
mon during the feedback loop. We address this point 
quantitatively below. 

An ideal reset function always prepares the same pure 
qubit state, regardless of its input. To fully quantify the 
performance of our reset scheme, we measure its effect 
on different superposition states \9). After applying Fb 
(Fbi), we extract P| ) by averaging over 300 000 repe- 
titions of the experiment (Fig. 3). Without feedback, 
P|o) is trivially a sinusoidal function of 9, with contrast 
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FIG. 3. (color online). Deterministic reset from any qubit 
state. Ground-state population P| ) as a function of the initial 
state \9) = cos(#/2) |0}+sin(#/2) |1), prepared by coherent ro- 
tation after initialization in |0), as in Fig. 2. The cases shown 
are: no feedback (circles), one instance of Fbo (squares) or 
Fbi (diamonds), and two instances of Fbo (upward triangles) 
or Fbo followed by Fbi (downward triangles). Inset: error 
probabilities for two rounds of feedback, defined as 1 — P\ t ) , 
where \t) £ {0, 1} is the target state. The vertical axis is 
calibrated with the reference states (Fig. 2) and corrected for 
the partial excitation between Mo and Mi. The systematic 
~ 0.3% difference between the two cases is attributed to error 
in the n pulse preceding the measurement of P|i) following 
Fbi [22]. Curves: model including readout errors and equili- 
bration [22]. 

slightly reduced from unity due to a preparation infi- 
delity of 0.7%. With feedback, we observe a suppression 
of the Rabi oscillation, with P| ) approaching the ideal 
value 1 (0) for Fb (Fbi) for any input state. However, 
a dependence on 9 remains, with P err = 1 — P| ) for Fbo 
(1 - P {1) for Fbi) ranging from 1.2% (1.4%) for 9 = 
to 7.8% (8.4%) for 9 = ir. The remaining errors have 
two sources: mismatch between measurement result M 
and post-measurement state occurring with probabil- 
ity pf* , for initial state and equilibration during the 
t = 2.4 fis lapse between the end of Mi and the start of 
M2, set by processing time in the controller. Transitions 
to 1 2) during Mi (with probability P21 = P21 + P2i)i or 
during r (r2ir), cause leakage out of the qubit subspace, 
where the feedback has no action. For perfect pulses, the 
overall errors (equal for Fbo an d Fbi) are to first order: 

^rr°=P00+Pf0+ri0T, 

P^" = P11 + P01 +P21 + (Toi + r 2 i)r, 

and weighted combinations thereof for other 9. Using 
the best-fit rates F^ (Fig. 1), errors due to equilibration 
sum to 0.7% (6.9%) for 9 = (n), while readout errors 
account for the remaining 0.4% (1.4%). 

To improve reset fidelities, we concatenate two feed- 
back cycles. The dominant error for 9 = remains un- 
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FIG. 4. (color online). Fast qubit reset. Initialization errors 
under looped execution of a simple experiment consisting of 
(a) a 7r pulse and (b) no pulses, leaving the qubit ideally in 
|1) or |0), respectively. Subsequent runs are separated by 
an initialization time Ti n j t . Black: initialization by waiting 
(no feedback). Purple: initialization by feedback, with three 
rounds of Fbo and a ir pulse on the 1 2 transition [22]. 
Curves correspond to a master equation simulation assuming 
perfect pulses. 



changed, but for 9 = tt it decreases to P cr 7° +P21 + ^2\t, 
in agreement with the measured values of 1.3% and 3.4%, 
respectively. The consecutive application of Fbo and Fbi 
prepares |1) with similar fidelities. 

The key advantage of reset by feedback is the ability to 
ready a qubit for further computation fast compared to 
coherence times available in 3D cQED [13, 28]. This will 
be important for ancilla qubits used for multiple rounds 
of error correction [29]. To test the gain in speed over 
passive initialization, we run a sequence multiple times, 
separated by a variable initialization time Tj n i t , during 
which we perform reset or not (Fig. 4). A measurement 
pulse follows the initialization period to quantify the ini- 
tialization error P err . The sequence itself emulates an al- 
gorithm leaving the qubit in |1) (or |0)), by simply apply- 
ing (or not applying) a tt pulse. Without feedback, when 
the sequence is a n pulse [Fig. 4(a)] the qubit remains 
partially excited and P orr approaches 50% as Tj n j t — > 0. 
With no 7T pulse [Fig. 4(b)], P err simply equals the total 
steady-state excitation. Feedback during rj n i t suppresses 
the reset error from ~ 32% (averaged between the two 
cases) to 3.5% for the shortest Ti n i t (15 /xs). We achieve 
this improvement by combining three rounds of Fbo with 
a 7T pulse on the 1 O 2 transition before the final Fbo, 
which partially counters leakage to the second excited 
state [22]. In the opposite limit Tj n i t ^> Ti, residual ex- 
citations are reduced from 16% to 3% [30], cooling the 
transmon. 

In conclusion, we have demonstrated feedback control 
of a transmon qubit using high-fidelity projective mea- 
surement and conditional operation. We have applied 
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this feedback to deterministically reset the qubit, starting 
from any superposition, with an average error of 2.4%. 
We have also used feedback-based initialization to loop a 
test sequence with more than 15-fold increase in repeti- 
tion rate relative to the passive method. Reset fidelity is 
currently limited by transitions to higher energy levels, 
caused by the 127 mK transmon temperature. We esti- 
mate that by reducing this temperature below 50 mK, 
the suppression of Tio and T 2 i will decrease the reset 
error past the ~ 1% fault-tolerant threshold in surface 
coding [31]. This reduction may be achieved by a combi- 
nation of infrared radiation shielding [24], use of a cop- 
per cavity [32], and improved qubit thermal anchoring. 
Moreover, decreasing the latency in the feedback loop, 
for example by using field-programmable gate arrays, will 
further reduce errors due to equilibration between mea- 
surement and action. While this demonstration employs 
feedback on a single qubit, the scheme can be extended 
to conditionally drive another qubit within its coherence 
time, realizing the feedforward [1] needed for teleporta- 
tion and measurement-based error correction. Future ex- 
periments will also target the generation of entanglement 
by combining the demonstrated feedback scheme with 
cavity- based parity measurements [33, 34]. 
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R. N. Schouten, and A. Wallraff for discussions. We 
acknowledge funding from the Dutch Organization for 
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Note added. — A parallel manuscript from Ecole Nor- 
male Superieure, Paris [35] reports feedback control of a 
transmon qubit using a Josephson parametric converter 
for high-fidelity projective readout. 
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READOUT AND FEEDBACK 
CONFIGURATIONS 

The readout frequency is set to wrf = w r — \ to maxi- 
mize the phase difference in cavity transmission between 
qubit in |0) and |1). The JPA is tuned to wrf by the 
combination of external flux biasing and a continuous- 
wave pump tone that bends the JPA resonance, as de- 
scribed in Rcf. 1. In contrast to this previous work, here 
we use a single generator for pump, measurement and 
local oscillator in order to maximize phase stability. A 
phase shifter on the RF arm is used to maximize the 
phase-sensitive amplification [Fig. S3(a)]. Another phase 
shifter at the output side maximizes the in-phase signal 
after IQ demodulation. An additional circulator between 
cavity and JPA is used to further suppress leakage of the 
pump into the cavity (Fig. SI). 

The ADwin samples the readout signal once, at a vari- 
able delay following a trigger from an arbitrary waveform 
generator (Tektronix AWG5014). This delay is adjusted 
in 100 ns steps to maximize readout contrast between 
states prepared with and without a 7r pulse on the Ool 
transition. For multiple feedback instances, the timing 
for each measurement is separately optimized. Once 
these timings are set, a routine (repeated every ~ 5 min), 
determines the optimum threshold for digitizing the read- 
out signal. This voltage is then used to assign H or L to 
the measurement in any feedback cycle. For Fb (Fbi), 
the ADwin triggers another arbitrary waveform genera- 
tor (Tektronix AWG520) to produce a it pulse when the 
outcome is L (H). Pulse timings and signal delays in the 
feedback cycle are illustrated in Fig. S2. 

MEASUREMENT OF THE TRANSMON 
POPULATIONS 

Whereas readout for feedback maximally discriminates 
between |0) and |1) from single measurements, the final 
readout, used to quantify the reset error P e „, relies on 
the average of thousands of runs. Due to leakage from 
the qubit subspace, the output voltage Vh also depends 
on P| 2 ) . The average measurement is related to the pop- 
ulations as [3] 

(V H > = (V[o>) P| > + <^i>) P|i> + (V m ) P|2>, (SI) 



where the (Vjj)) are the mean values of the distributions 
in Fig. S3. When (Vji>) ^ (Vj 2 >), as in Fig. S3(a), the 
measurement of (Vh) is not sufficient to extract the three 
populations. Two possible strategies to obtain two ad- 
ditional equations are: 1) repeating the sequence, but 
swapping two of the populations prior to measurement in 
order to obtain different (V^), (V{[); 2) analyzing the full 
time-domain response of the cavity [3] . We use a different 
approach, taking advantage of the degree of freedom in 
the relative phase between RF and pump, to extract P| ) 
with the minimum number of measurements. As shown 
in Fig. S3(b), the phase between RF and pump is chosen 
to maximize overlap between the histograms for |1) and 
|2). Making (Vji>) = (Vj 2) ), and assuming ELo P K> = 1 
(as justified below), simplifies Eq. (SI) to 

(Vh)=P|o><%) + (1-P|o>)<V| 1) ), (S2) 

with P| > as the only unknown. To measure Pji), we 
transfer this population to |0) with a it pulse, making 
(Va) = P\i) (Vfo>> + (f - Pji>) <Vfi>>. Similarly, to mea- 
sure P| 2 ), we transfer this population to |0) with two con- 
secutive 7r pulses, the first on 1 -H- 2 and the second on 
Ool. The amplitudes (Vjo)) , (Vji)) and (V\2)) are cali- 
brated by preparing the corresponding states. To do this, 
we initialize to |0) as best we can (P\o) — 99.3%) by per- 
forming a measurement (used for postsclcction [1, 4]) and 
then applying the required n pulses. Values for (Vj )), 
(Vm) and (V| 2 )) are estimated from <~ 10 000 shots each. 
As shown in Fig. S3(b), (Vji)) and (Vj 2 )) are matched to 
better than 1%. The sum of the extracted populations 
shown in Figs. S4(a-b) averages to 99.7% and 99.5% over 
the temporal evolutions starting from |0) and |1), re- 
spectively. For initial state |2) [Fig. S4(c)], leakage to 
the third excited state is relevant under 50 /is. In most 
cases, however, P| 2 ) is limited to a few percent, justifying 
the truncation of the model to the lowest three transmon 
levels. 



MODEL FOR MEASUREMENT AND 
FEEDBACK 

The readout errors ef 1 , indicating the probability of a 
wrong result M for the pre- measurement state \i), are 
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FIG. SI. Experimental setup. The key components of the feedback loop are highlighted in red. A Tektronix AWG5014 generates 
modulation envelopes for measurement and deterministic qubit pulses. The readout signal is amplified with a JPA at base 
temperature, a HEMT (Caltech CITCRYO 1-12) at 3 K, and two GaAs FET amplifiers (Miteq AFS3) at room temperature 
prior to demodulation to baseband (0 Hz intermediate frequency). The demodulated signal is amplified (Stanford Research 
Systems SR445A) and then split. One arm is sent through a low-pass filter (Mini-Circuits BLP-5+, DC-5 MHz passband) to an 
analog input of an ADwin Gold processor, and the other to the data acquisition card ( AlazarTech) . The ADwin compares the 
input voltage to the optimum threshold for discriminating between |0) and |1). Depending on the feedback logic programmed 
into it, the ADwin sends a trigger to the Tektronix AWG520 when the input is above/below the threshold. The triggered 
AWG520 produces the modulation envelope for a it pulse sent via the input line, closing the feedback loop. 
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FIG. S2. Feedback loop latency. The AWG5014 modulates the RF tone to generate a 400 ns measurement pulse reaching 
the cavity at t = 0. The black curve is a sketch of the photon population inside the cavity, which reaches ~ 45 photons at 
t = 0.4 /xs. We estimate a critical photon number n cr it w 150 and inverse measurement rate T" 1 = 26 ns at steady state [2]. 
The ADwin, triggered by the AWG5014, measures one channel of the output homodyne signal (red: qubit in |0), blue: |1)), 
delayed by ~ 200 ns due to the low pass filter at its input side. After comparison of the measured voltage at t — 0.6 /is to the 
reference threshold, the AWG520 is conditionally triggered at t — 2.54 fis, resulting in a n pulse (Gaussian, a = 6 ns, total 
length 24 ns), reaching the cavity at 2.62 /is. The ADwin acquisition (and the triggered -k pulse) can be shifted by ±0.1 /is to 
maximize readout fidelity. The measurement pulse for another feedback cycle or final readout reaches the cavity at t — 2.8 /is. 
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FIG. S3. The two readout configurations employed and the 
associated histograms for the three lowest-energy transmon 
levels, (a) Sketch in phase space of the coherent states cor- 
responding to the first three transmon levels. The amplified 
quadrature, chosen by tuning the RF-pump phase difference, 
is either I or I', giving the histograms in (b) and (c), re- 
spectively, (b) Optimum discrimination between |0) and 
This choice is used in Mo and in all the measurements which 
are part of a feedback sequence to maximize the single-shot 
readout fidelity, (b) Readout configuration making |1) and 
1 2} indistinguishable, as their respective coherent states have 
equal projection onto the I' axis. This configuration is used in 
the final measurement of every sequence to accurately extract 
the level populations, as described in the text. 



extracted from the histograms in Figs. 2 and S3, using 
the result of the fit in Fig. S4 to account for the transmon 
equilibration between M and Mi [1] . The data yield the 
values e£ = 0.3%, ef = 0.7%. By analogy with Ref. I, 



P01 



we approximate the errors with e{{ r 
(notation defined in the main text) 
with the best-fit values for Fig. S3, 
and + Pqi = 1.4%, we estimate the remaining errors 



Poo and e i 
Combining these 

Poo + P10 = 0.4% 



pf = 0.1%, p L m = 1.0%. 



COMPARISON OF FEEDBACK PROTOCOLS 

A double round of feedback tends to equalize the ini- 
tialization error for different input states (Figs. 3 and S5). 
Augmenting the protocol to three or more rounds only 
marginally affects the residual difference, which is mainly 
due to leakage into |2). When a n pulse is repeatedly ap- 
plied on the 0f>l transition [Figs. 3(a) and S5(a)], this 
leakage depends on the fraction of time that |1) is pop- 
ulated, which increases at shorter Ti n i t . Therefore, P\2) 
can exceed its steady-state value, leading to increased 
P err . We address this problem by inverting the popu- 
lations in |1) and |2) after two iterations of Fbo- This 
brings the population, previously in |2), within reach of 
a third instance of Fbo- This strategy is only effective 
when P| 2 ) > P\\) , explaining why a residual error of ~ 3% 
is observed at long ri n ; t , when, after two rounds of feed- 
back, P|!) ss P\2)- Similarly, the repeated initialization 
in |0) without periodically applying a tt pulse [Figs. 4(b) 
and S5(b)] limits the escape probability to the excited 
states, as evident from the reduced error at short t^h- 
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FIG. S4. Transmon population dynamics showing residual excitation. Populations (a) P\o), (b) P\i), and (c) P\2) as a function 
of time following preparation in |0) (red), |1) (blue), and |2) (green). The combined fit of the master equation model to these 
data gives the transition rates and steady-state populations shown in Fig. 1. 
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FIG. S5. Fast reset using multiple instances of Fbo. The 
measurement sequence is illustrated in Fig. 4, with (a) 6 — tt, 
and (b) 9 = 0. At any Tinit Pen decreases by a factor ~ 3 
with one instance of Fbo compared to passive initialization 
(Fig. 4). Two or three iterations add only a modest improve- 
ment. An additional decrease in error by a factor ~ 2 at 
short ri n it is obtained with a n pulse on the 1 O 2 transition 
preceding the third instance (purple). This pulse corrects for 
the leakage of population to |2), occurring for short That. The 
last feedback round brings this population back to |0). Solid 
curves are given by the model, using the parameters extracted 
from Figs. 1 and 2 and assuming perfect pulses. 



